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The response of an open spin cliain witfi isotropic antiferromagnetic interactions to a uniform mag- 
netic field is studied by classical Monte Carlo simulations. It is observed how the induced on-site 
magnetization is non uniform, due to the occurrence of edge staggered terms which decay ex- 
ponentially over a distance equal to the zero field correlation length of the infinite chain. The 
total magnetic moment associated to each staggered term is found to be about half of the original 
single-spin magnitude and to decrease as the inverse of temperature (i.e. to behave as a Curie-like 
moment). The numerical results are compared to recent NMR findings in spinless-doped Y2BaNi05; 
the remarkable agreement found shows that, for temperatures above the Haldane gap, the classical 
approach gives a correct picture of the boundary effects observed in the Heisenberg S=l chain. 

PACS numbers; 75.10.Jm, 75.40. Mg, 75.50.Ee, 76.60.-k 

Because of the richness of their phase diagrams, 
Heisenberg spin systems with low-dimensional antifer- 
romagnetic interactions (HAF's) are presently attract- 
ing strong interest. Magnetic correlations in spatially- 
homogeneus HAF's have been probed dynamically, in 
Fourier-space, by neutron scattering, but can hardly be 
visualized experimentally in a static real-space picture. 
Translational invariance, infact, makes each site equiva- 
lent to the others, thus preventing any spatial oscillation 
of the spin direction. In recent experiments, the cor- 
relation properties of low-dimensional HAF's have been 
investigated through NMR imaging of the spin polariza- 
tion induced by a uniform field around non-magnetic de- 
fects that break the translational invariance. Relevant re- 
sults have been obtained in cogper oxide two-dimensional 
compoundsElu, in spin laddersH, in half-integeiu and inte- 
ger spin chainsEl. 

In S" = 1 chains in particular, the local magnetization 
{Sf) has been resolved psife hy site, through NMR 
in Mg-doped Y2BaNi05la, for temperatures ranging from 
T ~ 0.35 J/ A:b to T ~ J/fcs, being J the AF exchange 
constant. (S'f ) shows an alternate component which is 
maximum around impurities (i.e. close to chain bound- 
aries) and vanishes exponentially over a distance equal to 
the zero-field correlation length of the bulk. Boundary 
staggered defects with total spin S — 1/2 and a size of the 
order of the bulk correlation length are actually expected, 
ioT S — 1 HAF chains, in the limit T 0. At very low 
temperatures in fact, the magnetic properties of t 
systems are controlled by edge- induced triplet, statesi 
in which {Sf) shows the profile described aboveBlB. This 
T ~ argument however, can hardly be used by itself 
in order to explain the experimental evidences in Ref. 
5, since - on increasing temperature - the excited states 
above the Haldane gap A/f ~ 0.4J should be taken into 
account (in open S — \ chains A^f corresponds to the 
separation between the second and the third lowest ly- 



using Quan- 



which all the excited states are corrects 
recently carried out by Alet and S0rensenc. 
tum Monte Carlo techniques. Excellent agreementrjwith 
the experimental data derived from NMR imagingQ has 
been obtained. 

The purpose of the present work is to study whether 
or not the alternating magnetization observed in finite 
S — 1 chains is a pure quantum mechanical effect some- 
what reminiscent of the liquid-like ground state. The re- 
sponse of open chains with various lengths L to a uniform 
magnetic field is thus studied in the limit of infinite-S*, 
proving that edge staggered defects with spin-S/2 origi- 
nate from the translational invariance breaking, even in 
the framework of a classical model. It is shown that 
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FIG. 1. Magnetization profiles at T = 150 K, obtained 
by classical Monte Carlo in a 5* = 1 chain of 31 sites with 
J/kB~2S5 K. Different symbols correspond to different ap- 
plied fields: H = Q (stars), H = 5 Tesla (triangles) and 
H = 14.1 Tesla (squares). 
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FIG. 2. Left side: magnetization profiles (Sf) in a field H — 14.1 Tesla for 5 = 1 chains with J/fcs ~ 285K and various 
lengths (from classical Monte Carlo). Solid lines are fittings of the data according to the phenomenological model described by 
Eq. @. On the right side NMR spectra in doped HAF chain Y2BaNii_^Mg^05 (from Ref. 5) are reported. Through Eq. 
(^) in the text, the local magnetization of the classical chain can be quantitatively compared to the NMR spectra of the S—l 
chain and the correspondence between boundary moments and satellite peaks is directly inferred. 



the characteristic spatial extention of these defects corre- 
sponds to the zero-field correlation length of the thermo- 
dynamic chain {L — s- oo) and that, at high temperature, 
the behavior of (5f ) calculated for infinite-S' tracks the 
experimental findings in the S — I chain. 

In the classical limit (S* — + oo), the Hamiltonian for a 
nearest-neighbor Heisenberg chain of L sites in an exter- 
nal field H takes the form 



(1) 



where S'i's are classical vectors, whose magnitude is taken 

\/2 in this case), and 



to be 



S^ 



y/SjsTT) (i.e 



J is positive for antiferromagnets. 

The thermal expectation value of the on-site magnetic 
moment has been calculated by the Wolfili3 cluster algo- 
rithm in the temperature range 5QK < T < 285K, for 
L up to 31 spins. In order to compare numerical data 
with NMR results in. doped Y2BaNi05, J/ks has been 
set equal to 285 Kllil. The effect of the magnetic field 
H has been included using a standard Metropolis algo- 
rithms, which associates a flip probability to the whole 
cluster. As the field intensity is low, the flip acceptance 
remains close to one. 

Typical results of the simulations for different values 
of H are shown in Fig. |l|, proving that the response of a 
finite classical chain to a homogeneous field is non uni- 
form. In Fig. |2[ magnetization profiles at characteristic 



temperatures and i/=14.1 Tesla, are shown together with 
^^Y NMR spectra in Mg-doped Y2BaNi05 (at the same 
temperatures) from Ref. 5. The spectra were obtained at 
fixed frequency {I'rj =29.4 MHz) by sweeping the mag- 
netic field Hq in a narrow range around 14.1 Tesla. The 
intensity of the NMR signal is proportional to the num- 
ber of ^^Y nuclei which obey the resonance condition 



Hn = 27r 



- k 



A{St 



(2) 



where 7 is the ®^Y-gyromagnetic factor, A =1.3 Tesla is 
the *^Y - Ni^+ hyperfine coupling constant in Y2BaNi05 
and k is an i-independent factor that accounts for the 
chemical and orbital shifts. Using Eq. (||), the local 
magnetization of the open classical chain and the posi- 
tion of the peaks in the NMR spectra can be directly 
compared. Since k is not known precisely a priori, we 
fix it by matching the spin-polarization at the centre_af 
the chain with the maximum of the NMR central lineEj. 
Then, as sketched in Fig. ^ a satellite peak is found in 
correspondence of each value of magnetization taken by 
the edge spins. 

The magnetization profiles can be analyzed in detail, 
by observing that {Sf) consists in a uniform part and of 
staggered contributions decaying away from each bound- 
ary. The following phenomenological function 
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= {S^) + {S!) [(-l)(-i)e^ + (-l)(^-^)e^ 



(3) 



is thus used to fit the numerical data. Expression (|^) 
describes very well the behavior of (5*?) in all the in- 
vestigated temperature range (see for instance solid lines 
in Fig. |2|) and the extracted fitting parameters - (5*^), 
(Sf) and ^ - turn out to be independent of the chain 
length L. The two alternating terms in Eq. (^) give an 
in phase (out of phase) contribution in chains with odd 
(even) number of spins. In real systems, consisting in 
an ensamble of segments with different length, such ef- 
fect induces a distribution of (S^) and thus a broadening 
of the NMR satellite peaks, as actually observed in Fig. 

Here we do not investigate this subject quantitatively, 
being mainly interested to discuss the behavior of the pa- 
rameters extracted from the classical Monte Carlo data 
in light of the values for the same quantities obtained by 
NMR in the S^l chain. 

In Fig. ^), the temperature behavior of the uni- 
form term {S§) is reported. {S§) depends linearly on 
the field up to 40 Tesla (inset) and follows strictly the 
magnetization x ' H that would be obtained by using 
for X the_pero-field susceptibility of the infinite classi- 
cal chairll3 (solid line in Fig. ||a)). This result should 
be intuitively expected, since the staggered terms in Eq. 
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FIG. 3. a) Temperature beliavior of the uniform magneti- 
zation {S^} in H =14.1 Tesla extracted from classical Monte 
Carlo data by Eq. (^ (full squares). Solid line is the exact 
classical thermal expectation value of the on-site magnetic 
moment for a chain with L oo, calculated through Fisher's 
formula for the uniform susceptibilityliS. b) The first-site 
staggered magnetization (Si) of the classical chain is com- 
pared to the experimental value of (Sf) measured by NMR in 
Mg-doped Y2BaNi05Q (open circles). In the insets it is shown 
how (S^) and {Sf) depend linearly on the applied magnetic 
field up to 40 Tesla. 
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FIG. 4. Decay length ^ of the staggered magnetization, 
obtained from Monte Carlo calculations through Eq. (^) 
(full squaresV, vis a vis to the experimental values of ^ 
(open circles )u and to ths zero field correlation length of the 
infinite-L classical chainu (solid line). In inset, the total mag- 
netization (Sj-), associated to one of the chain boundaries (see 
Eq. (0) in the text), is fitted by the Curie-law (S). 



(^) involve only boundary spins and thus, in the ther- 
modynamic limit, the "bulk" magnetization is due only 
to (S'b). The first-site magnitude (Sf) and the decay- 
length ^ of the staggered contributions, extracted from 
the simulation, are reported respectively in Fig. jS^j) and 
^ and compared with the experimental v^alues of (S'f) 
and ^ from ^^Y NMR in doped YzBaNiOgQ. The agree- 
ment between experimental results in the HAF S~l one- 
dimensional compound and predictions of the classical 
model is remarkable at high temperature. Small quan- 
titative deviations, which likely prelude to a more sensi- 
ble departure at lower temperatures, are observed only 
around T =100K~ An/kB- At the moment, the lack 
of experimental data for (5J) and ^ below 100 K (due 
to the broadening of the NMR lines) prevents a precise 
definition of the temperature range in which the classical 
model accounts for the magnetic response of the open 
S' = 1 chain. Fig. ^ also shows how the decay length ^ 
of the staggered magnetization in a weak external field 
{H—lA.l Tesla) tracks strictly the behavior of the zero- 
field spiUpSpin correlation length (solid line), calculated 
by FisherliS for the infinite-volume chain. 

In the inset of Fig. ^, the total magnetic moment as- 
sociated to each staggered contribution. 



(4) 



is plotted as a function of temperature. (S^) is well re- 
produced by a Curie-like law 



(5) 



characteristic of non-interacting classical moments of 
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magnitude Sf. The value of Sf that optimizes the fit- 
ting (sohd Une in the inset in Fig. ^ is 0.72, about half 
of the site-spin magnitude ^/ S {S + 1) = \f2. This re- 
sult extends to a temperature region in which the system 
shows classical behavior the low-T picture for gapped 
spin chains, in which. |iSy2 degrees of freedom develop 
at the chain edgesEfra^EJ. Moreover, since the temper- 
ature dependence of is mainly controlled by 
for Ajj/fes < T < J/ks, even the first-site staggered 
moment displays approximatively a Curie-like behavior. 

In summary, our numerical results prove that the al- 
ternate boundary magnetization observed experimentally 
in finite-length 5=1 HAF chains is correctly predicted, 
for An/ks < T < J/ks, by a strictly classical anal- 
ysis of the Heisenbeg model. In light of this evidence 
we conclude that the occurrence of Curie-like fractional- 
spin defects is a general feature in HAF one-dimensional 
systems with open boundaries, not related to the occur- 
rence of a spin gap. On the other hand, the size of these 
staggered defects, which is fixed by the bulk spin-spin 
correlation length, is affected by quantum fluctuations 
and is thus expected to remain finite when the ground 
state is spin-liquid like. 

The authors are grateful to A. Rigamonti, E. S0rensen 
and M. Fabrizio for stimulating discussions. 
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Assuming 7=27r-2.0859 MHz/Tesla, the shift of the NMR 
central line due to chemical and orbital contributions, ex- 
tracted from the comparison between numerical and exper- 
imental results at T = 283 K, turns out k/vRp =390ppm, 
in good asMpment with previous extimations in the pure 
compouncO'Lj. 
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